Microstructure of He II in the presence of boundaries 



(N 

o 

(N 

> 
O 

o 



5-H 

CD 



i 

O 
o 



> 
in 

00 



o 

(N 



X 



M. D. TomchenkcB 

Bogolyubov Institute for Theoretical Physics, l^-b Metrolohichna Street, Kiev 03680, Ukraine 

We show that, in spite of the commonly accepted expectations, the boundaries affect strongly 
the bulk microstructure of He II. We determine the structure of exact wave functions of the ground 
(^o) and excited (*l/k) states of He II under zero boundary conditions. The explicit account for 
boundaries leads to a significant change of the equations for and and the theoretical values 
of the energy of quasiparticles E(k) and the energy of the ground state Eo, as compared with the 
solutions for periodic boundary conditions. In the limit of a weak interaction, we obtain the formula 



+ WfififM ( d i s the space dimensionality), which differs from the well-known 



Bogolyubov formula by the factor 2~ d . The effect of boundaries is topological and is related to the 
different collections of eigenmodes for closed and open system. 

PACS numbers: 67.10.-j, 67.25. dj, 67.25. dt 



I. INRODUCTION 

Starting from the pioneer works by L. Landau N. Bo- 
golyubov 0, [j| , and R. Feynman [3, Q , the microscopic 
theory of superfluid He A (He II) is developed for about 
seven decades. It explains sufficiently well the basic prop- 
erties of He II (see reviews |6H8J ) and is considered to be 
established. At the same time, some questions remain 
open. In particular, the nature of the A-transition, num- 
ber of higher condensates, properties and role of micro- 
scopic vortex rings are not clear. Moreover, the relation- 
ship of superfluidity to the presence of a condensate of 
atoms is not completely clarified. 

In addition, the construction of a microscopic model 
uses always the periodic boundary conditions (BCs) with 
the subsequent transition to the thermodynamic limit 
(N,V ~ > oo, N/V — const) instead of realistic BCs (a 
finite system with the surface). Certain difficulties in the 
construction of a perturbation theory at realistic BCs 
were indicated in the articles by A. Bijl Q and M. Gi- 
rardeau [l3| and the monograph by K. Huang This 
is due to the fact that the bare wave function (WF) is nat- 
urally sought in the form of a solution for the ideal gas, as 
a product of sines sin (jrXj/L x ) (in the one-dimensional 
case) . But such sines set a very inhomogeneous "arc-like" 
density of a medium, whereas the density of real helium is 
homogeneous. In view of this difficulty, the realistic BCs 
were changed by the periodic ones in order to simplify 
the calculations. 

To our knowledge, the solutions for the WFs of inter- 
acting Bose particles with the explicit account for bound- 
aries were obtained only for a one-dimensional system of 
particles with point interaction [Hj]. We will find 
the solutions for particles with any nonpoint interatomic 
potential possessing the Fourier-transform and for any 
dimensionality. 



It is clear intuitively that the direct influence of bound- 
aries on a bulk Bose fluid (e.g., He II) should be small. 
But He II is a strongly nonlinear and highly correlated 
system with a condensate. Since the condensate corre- 
sponds to a delicate tuning of the whole system, even 
a very weak effect of boundaries can change it. We will 
show that the boundaries strongly affect a quantum fluid: 
they affect the spectrum of eigenmodes of the ground and 
excited states. As a result, the calculated values of the 
ground state (GS) energy, energy of quasiparticles, and a 
structure of a condensate [14| become essentially different 
from those for periodic BCs. 



II. GROUND STATE OF HE II: EQUATIONS 

Several methods of determination of the WF of 
the GS of helium-II [I. ll5rl2(| (see also reviews 0,11]) are 
available. They used periodic BCs, which cannot be obvi- 
ously realized in the nature. Our task is to find under 
BCs close to those realized in the nature and for realistic 
interatomic potentials. The influence of boundaries on 
the structure of the WF can be nontrivial. Therefore, 
it is expedient to use one of the exact analytic methods 
for the determination of VPo [HI~Gil H3] j since the numer- 
ical methods may not catch the effect. We will operate 
within the method of collective variables [l?], H3| used by 
us [Hill. 

We transit to the collective variables: 



Pk = 



1 N 



(k^O), 



(1) 



where N is the total number of atoms, and Tj are the 
coordinates of atoms. For convenience, we introduce the 
following notation. The sums over the wave vectors 



2tt 



Jx Jy Jz 



(2) 
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L x Ly L 

{j x ,jy,jz are integers, and L Xl L y ,L z are the sizes of the 



2 



(27T) 

system) are denoted by the index (2ir); for example, ^ . 

k 

And the sums over the wave vectors 



Jx Jy Jz 
L X Ly L z 



(3) 



(tt) 



are denoted by the index (tt). For example, means 

k,q#0 

the sum over k and q, which are nonzero and run series 
([3]). If the row with (tt) or (2tt) contains some designation 
(e.g., k + q ^ 0), it has no relation to (tt) or (2tt). 

At the periodic BCs, the solution for is obtained 
by I. Yukhnovskii and I. Vakarchuk [17] : 



(4) 



x - a 2 (ki) 
bb = — 2! — PklP_kl 



(2-7r)ki+k 2 #0 



+ E 



a 3 (kx,k 2 ) 



ki,k 2 #0 
(27r)ki+k 2 +k 3 ^0 



3WN 



"PkiPksP-ki-ka 



(5) 



E 



ki,k 2 ,k 3 ^0 



a 4 (ki,k2,k3 
4!A 



" Pk! Pk 2 Pk 3 P- k i - k 2 - k 3 



+ 



where the functions Oj(ki, . . . ,kj_i) satisfy the chain of 
equations 



N-l i-^jWfc) h 2 k 2 „ . 
En = nv 3 (0) - > — — - > oofk), 6) 

k^O k#0 



2A 



nv^(k)m 



a 2 (k)fc 2 - a^(k)A; 2 = 



(7) 



(2tt) . . (2tt) 



q#0 



o 3 (k,q) w-2 



i?(k,q) 



e (fc) + e (9) +e (k + q) : 



(8) 



R(k, q) = kqa 2 (k)a 2 (q) - k(k + q)a 2 (k)a 2 (k + q) - 
- q(k + q)a 2 (q)a 2 (k + q), (9) 



e (k) = fc 2 (l-2a 2 (k)), 



(10) 



where n — N/V, V is the volume of the system, Eq is 
the "bulk" energy of the GS per atom, and 



i*(k) = V3(k) = / U 3 (r)e- lkr dr 



(11) 



is the Fourier-transform of the interaction potential U%(r) 
of two -He 4 atoms. The equation for 03 is written in 
the zero approximation not involving the corrections 
04 and 05, and the equations for higher corrections 
a 3 >4(ki, . . . , kj_i) are not presented. 

Consider He II in a rectangular vessel L x x L y x L z 
in size. In order to solve the problem exactly, the solu- 
tion should be sought inside and outside of the vessel, 
and the continuity of the WF and all its derivatives on 
the boundary can be used as the BCs. However, it is 
too complicated. To simplify, let us consider imperme- 
able walls. Then the WF is nonzero inside the vessel, 
and outside of it and on the boundary is equal to zero. 
The real walls are not impermeable, but they create a 
high potential barrier with U s > 1000 K (see Eq. (fT4|)). 
At such a barrier, the solution for a single particle in 
the box differs from the solution under an infinitely high 
(impermeable) barrier arbitrarily slightly. 

The Hamiltonian of liquid helium bounded by the walls 
has the form 



H 



2m 



E A ^ 



1 



2E 



Ua(\ ri 



i\) + U 2 + U em , (12) 



U 2 = ^ [Uw(xj) + U W (L X - xj) + U w (yj) + 
3 

+ U W (L V - yj ) + U w (zj) + U W (L Z - Zj )] . (13) 

Here, U2 is the interaction potential of helium atoms with 
all walls. The potentials of the walls are similar by form 
29] to the interaction potential of two helium atoms. The 
free surface of liquid helium creates an energy barrier 
also. Therefore, we consider it as a wall. Below, we 
assume all barriers identical with simple potential 



U w (x) 



U s >0 x<0, 
x > 0. 



(14) 



Below, we take U s finite, and the transition U s —¥ 00 will 
be made in final formulas. 

Works [13, H3] used the non-Hermitian Bogolyubov- 
Zubarev Hamiltonian. Since the absence of the Hermi- 
tian property induces the questions (30}, we will start 
from the ordinary Hermitian Hamiltonian (|12p . 

As is known, a helium atom disturbs the distribution 
of charges in a wall. If a wall is metallic, the perturbation 
is equivalent to an exact mirror "image" of the atom 
in the wall. In the image, the nucleus and "electrons" 
have, respectively, negative and positive charges. If the 
wall is made of a dielectric, the charges of the image 
are diminished. The interaction with images is taken 
into account in U em . However, the consideration of the 
images strongly complicates the calculations. Therefore, 
we set U em = in what follows. 

We expand the potential in a Fourier series: 



U 3 (\r 1 -r 2 \) 



1 



2 d V ^ 

k 



ZAj(k)e 



ik(ri— r 2 ) 



(15) 



3 



where d = 3 is the space dimensionality. We seek the 
WF of helium atoms inside the vessel, since the WF 
is zero outside of it. In this case, all functions in the 
Schrodinger equation should be expanded in a Fourier 
series with regard for the fact that the coordinates Vj of 
atoms are defined only in the interior of the vessel. As is 
known, such an expansion reproduces exactly the values 
of functions for all inside of the vessel. In the poten- 
tial £/3(|rj — rj\), the coordinates Xi and Xj take values 
from the interval [0,LJ, whereas the values of %i — Xj 
belong to the interval [— L X ,L X ]. Analogously, we con- 
sider the y- and z-components. By the rules of Fourier 
analysis, the function Us(\ri —Tj\) must be expanded in 
this case in series (fTS"|): with k (J3j> and the volume equal 
to 2L X x 2L y x 2L Z = 2 d V. 

For the total interatomic potential, we obtain 



1^,, ,| ,x ra/ a (k) 



k#0 



+ (N-l) 



ni/ 3 (0) 
2 d+i 



E 

k#0 



2 d+i ■ 



(16) 



Under periodic BCs, xi — Xj take values from the interval 
[0, L x ]. Therefore, we need to replace 2 d V in (fT5|) by V 
and to quantize k by law ([2]). Then we have Q 



I (2t) ,,s 



2 



k#0 ^ 



(17) 



This is a key point. Since Eqs. ©-(|H]) are "generated" by 
the total potential, we may assert at once that the explicit 
account for boundaries will lead to the replacement (in 
these equations) of Us (k) by v^{k)/2 d and the summation 
law © by ©. 

Let us determine the WF for the system of Bosc 
particles in a vessel. vPq satisfies the Schrodinger equation 



(18) 



with Hamiltonian (jl^ - flTBl . where we neglect U em . 

The main question is in which form we should seek 
* - The solution of gj), © for * with periodic BCs is 
usually obtained with regard for the requirement of the 
translation invariance and the functional independence 
[l?} of the following collections of variables p^: 



P-ku P-kiP-ka 



P-kiP-k 2 P-k 3 



(19) 



But the boundaries break the translation invariance. 
However, the solution for ©, ([5]) can be obtained with- 
out regard for the translation invariance, by substituting 
the bare function — const, describing the free bosons, 
in the Schrodinger equation. The equation itself "sets" 



the form of corrections, which appear in the exponent. 
For the problem with boundaries, we will seek anal- 
ogously, by substituting the solution for free particles in 
the Schrodinger equation. 

The WFs of the system of N free Bose particles located 
in a rectangular box with potential (|13p , (|14[) are known. 
If the states of all particles are identical, the WF is equal 
to one of the functions [3l| 

N 

*f = n{ sin ( fc ^' + s ix) shx ( k i y yi +K) ahx ( k h z j + <U}> 

3=1 

(20) 

where l x ,l y ,l z — 1,2,3,..., and fc; x and Si x satisfy the 
equations 

. p Jxki x L x Jxki x L x 
sm ^x = i , sm(ki x L x + 5iJ = ± 



2 ' 
(21) 

which yields 

ki m L x = nl x - 25i x , Si x = arcsin {^ x ki x L x /2). (22) 



Here, we denote j x = YT\j an d the values of arcsin 

are taken between and 7r/2. For the realistic systems, 
7 X -C 1. For example, for He A atoms at U s = 100 if and 
L x = 10 cm we have ^ x ps 4.92 x 10~ 10 . At l x -C I/t^Yx, 
relation (|2"2"|) yields 



h x fa l x kix = ^(1 - lx)^/L x 



5l x ~ l x hx = lxTrj x /2. 



(23) 



(24) 



The maximally possible values of <5; x and ki x are deter- 
mined from the condition sin(5; x = 1: 



°i x =7T /2, k, =k lx /d lx . 



(25) 



The WF of the ground state of a single atom inside of 
the box, 

*o,i( r ) = sin (faxX + S lx ) sin (k ly y + S ly ) sin (k lz z + S lz ), 

is positive and nonzero. Outside of the box, the WF 
decays exponentially. 

With regard for the above consideration, the WF of 
the ground state of N interacting Bose particles located 
in the box can be sought in the form 



* = Ax%& 



(26) 



where is the WF of the ground state of N free Bose 
particles in the box (function (|20[) with l x = l y = l z = 

1), e Sb is a WF of the form Q, © with the sums £ 

instead of (this is seen from expansion (fTB"])). and e «° 
describes the cross terms (from and e Sb ) arising at the 
substitution of ^Qe Sb in the Schrodinger equation. By 



4 



the theorem on nodes, if the GS is nondegenerate, then 
is nonzero in all points in the box. Therefore, it is 
convenient to lift the function Sw to the exponent. The 
factor e Sb appears due to the interatomic interaction, like 
that in (fJJ. We note that \1/q will ensure the fulfillment 
of BCs and the transition to the WF of free particles at 
the switch-off of the interaction. This gives the structure 

of (251. 

Notice that the structure (|26|) was proposed previously 
by S. Yushchenko [S3]. 

The function Wo (|26|) must satisfy the Schrodinger 
equation (TIBf and be zero at the boundaries. The last 
property is ensured at U s — > oo by the factor Wq which 
satisfies the equations 



— £ A^ f + U 2 ¥ = NE^l (27) 



Whence Ci x (0) = 0, and the function Ci x (j x ) is approx- 
imated by the formula 



C lx (j K >0)«-i+^^ln| 1 



4:j X 5l 



(34) 



which is valid at j x < j m = l/5\ x with an error of < 3%. 

We note that cot (k\ x x 4- ^i^) satisfies the require- 
ments to functions expanded in the Fourier series. At 
Six = 0, we obtain cot (irx/ L x ), which does not satisfy 
these requirements, since the integral of the modulus of 
this function over the interval [0, L x ] is equal to infinity. 
In view of this difficulty, we preserve all S to be nonzero 
in calculations and then set S — > in final formulas. 

Using the difference of the functional structure of the 
collections pk (fT9]) . we determine the structure of the 
solution for from flU])-© and ©: 



E, 



h 2 h 2 

/ 1 1,2 _ ,„2 , ,„2 , ;„2 



2m 



') fcj. — k lx + k ly + k lz 



(28) 



The function Sb satisfies the equation 



2m 



'£A j e s » + -'£u 3 (\r i -T j \)e s » = NE b e s \ (29) 



Then Eq. (fl"8]) is reduced to the equation for S 



(i). 



2m ^ 



/ 



A^W + (V,5W) 2 + 2V^W^o + 



(30) 



= N(E -E b -E£). 
We note that 

~~jj"T a = ix^lx COt (fcia-Xj- + 5ia.) 



(31) 



+ iykiy cot (fciyj/j + + i z fci z cot (k\ z Zj + 8i x ). 

By seeking Wo inside of the vessel, we set cot {k lx x + 8\ x ) 
on the interval [0, L x \. In order to apply the method 
of collective variables, we expand the cotangent in the 
Fourier series: 



(2n) 



cot (k lx x + 5 lx ) = Ci x (j x )e l27V ^ x 



/ix 



(32) 



where j x runs all integers. We have 



Cl x {j x ) = —Ci x (—j x ) 



1/2 



dx 



cot (ki x x + S lx )e iq * 



o 



-2i J dx sin (2itj x x) cot (irx ~ 25\ x x + Si x ). (33) 




(27T) 



q+qi#o (i) 

q#0 q, qi #0 ViV 

q+ qi +q 2 #0 (l) 

s 3 (q,qi,q 2 ) 



q,qi,q 2 ^ 



where q and q^ run values ((2|) and ((3|), respectively (the 
first fact follows from (|3"2"|) . and the second does from 
(fT5j) and (J7J). The functions satisfy the symmetry 
relations 



S , W(q,q 1 ) = 4 1) (q,-q-q 1 ) 



5 , 3 1) (q>qi.q2) = 4 (q^qi 



(36) 



^ Hq,qi,-q-qi -q 2 ) (37) 



^ 1) (q,-q-qi-q 2 ,q 2 ), 



(i) 



and the analogous ones can be written for higher S 

The equations for the functions Sj^ can be found from 
Eq. (j30]) , if we substitute all functions in it and collect all 
terms referred to each collection in (Ti"9|) and the constant. 
In view of the functional difference of collections in (|19p , 
we equate these groups of terms to zero. This yields the 



chain of equations for Eq and the functions S 



E = E* + A 1 , 



(i). 



(38) 



At 



2m 



(2tt) 



n'-^E^W^^)- ( 39 ) 

) 

-j= ^k 1 xq x C lx {q x )S[ 1 \-ci x ) + {x -> y,z) 



(2tt) 



5 



^ 1} (q)eo(q) 



2 

N 



iV N2ki x q x a 2 (-q x )C lx (q x )8 cli q x + (x ->• y, z) + 

(tt) 

E {(?? + qiq)4 1) (q 5 qi)+ 



qi#o 



3g?4 1) (q>qi.-qi) + 
7V2< Zl 2 ^ 1) (q 1 )4 1) (q-qi,qi) + 

.5%/7V(^ - q 1 q)^ 1 )(q 1 )5( 1) (q- qi )} 

(2t) 



+ 



E 2ik Xx C\x(j>x) (q x -p x )Sf'(q- pj- 



^(q, q 1 )[e (q 1 ) + e (q + q 1 )] + 



(40) 



+ 2S< 1) (q)a 2 (-q 1 )qq 1 - g 2 Sf 3 (q)a 3 (q, q X ) = 

= <5 q ,q x \/^V'«fci 2 ;C'i 2; (g x ) {2qi :E a2(q 1 ) - ^03^, q x )} + 

(2tt) 

+ E ifcixAx(fc) 4(g x + g lx -^^(q- Px^i)- 



3 

(tt) 



6PxS ( 3 >(q-p x ,q 1 ,-q-q 1 ) + (x->y,z) 



1 f 2 

~7^f E 1 ^7(92 -qiq 2 )5 , 3 1) (q.qi -q2,q 2 )- 



q 2 (qi + q 2 + q)5 , 3 1) (q>qi J q2) 



- 4q 2 ( qi - q 2 + q)5p ) (q 2 )^ 1) (q - q 2 , qi ) + 
+ 6( Z2 2 ^ 1) (q 2 )5( 1) (q - q 2 , q 1; -q - qj + 
+ 4(q 1 +q 2 ) 2 5< 1) (q 2 ,q 1 )5* 1) (q-q 2 ,q 1 +q 2 )}. (41) 

For brevity, we do not give the equations for the higher 
functions 5j> 3 (q). Equation (|41[) should be symmetrized 
to satisfy For the validity of relation ([3"Uj) . it is suf- 
ficient to consider only Sj with the first even argument. 

We have the rule: the first argument of the functions Sj 
is always quantized by the even law (J5|)- Otherwise, we 
set sf ) = 0. 

Analogously, (j29l) yields the chain of equations for Eq 
and the functions aj : 

~h N-l . , i^rw 3 (ft) h 2 k 2 „ N , x 
^ = ___n^ 3 0)- > -r^-- > a 2 (k), (42 



nu 3 (k)m 
2 d H 2 



2 d + 1 N ^ 2mN 



a 2 (k)fc 2 -a 2 (k)fc 2 = A 2 (k) + (43) 



g^ (g2+kq) + g g2 Mq^) ; 



a 3 (k,q) PS A 3 (k,q) - 



2i?(k,q) 



eo(fe) + eo(g) +e (k + q) 



(44) 



These equations can be found from ([H])-© by the 

(2tt) (tt) 

changes 1/3 (ft) — > v 3 (k)/2 d , —> ^ and the addition 
of "surface" corrections Aj . The latter are related to the 
appearance of the terms 



2m 



E A2 ( k : 



PkP-k, 



k#o 



ft 2 



(7r)ki+k 2 #0 



-r- E ^3(ki,k 2 )p kl Pk 2 P-k 1 -k 2 , ■ • ■ 
2m ^— ' 

ki,k 2 #0 

in Eq. ([30]) . which have the structure of terms of Eq. 
In particular, 



(2*) . 



-4, 



iki x Ci x {-q x ) 



4(q, + ft,)4 1) (q :c ,k) 



+ 6qxS { 3 1] (q x ,-q x ,k) + (x^y,z) + 
1 (27r) r 

+ ^EK + qkisf'Ns^q.k)- 



4(k + q) 2 5( 1) (q,k)^(-q,-k) 
69 2 5| 1) (-q)5^ ) (q,k,-k)' 



(45) 



Though the majority of terms in ([29)) and ([30)) has the 
structure of the collections /Ok, which is inherent only in 
the given equation f ([2^|) or ([3D])), a small part of terms in 
(f3"U| has the structure of terms in (|29p . The criterion for 
the construction of the chains of equations for Sj and 
aj is the functional difference of the collections pk (|19[) . 
Basically, the structures of the terms in ([29]) and ([30]) 
are different. Therefore, the partition of Eq. (jT5J) into 
([29]) and ([30)) is justified and simplifies the determination 
of Sw ■ The collections pk, which are present in both 
equations, should be analyzed with the use of the initial 
equation (JT5J). We can verify that ^0 (j2"6")) satisfies Eq. 
(fT8[) . if aj and Aj satisfy Eqs. |42[) - ([4%[) . 

As U s ->• 00 and 5 -> 0, Eqs. ([23, O-gS]) are 

the required solution for the WF of the ground state of 
N interacting Bose particles filling a rectangular vessel. 

The analysis of Eqs. (|40[) and ([41]) indicates that the 
functions S^(q) and S2 (q, q x ) can be considered "one- 
dimensional" : 



s[ 1 \q)^S[ 1 \q)(6^ x 



J q,q, 



^2 X) (q>qi) ~ ^^^(q^qi 



a)/ 



q^O 



q^O 



<5q,q y (q y , qi ) + 8^ S^> (q, , Ql X47) 



+ <W) 



(1)/ 



(46) 



6 



(q = (q x ,q y ,q z ) = q^, + q y + q z ), since the nonone- function e Sb sets a homogeneous distribution of atoms, 
dimensional with q = (q x ,q Vl 0) (with possible per- 



mutations of the x-, y-, and z-components) and q = 
(q x , q y ,q z ) are less than one-dimensional ones by, respec- 
tively, iV 1 / 3 and TV 2 / 3 times. Due to the smallness of 
nonone-dimensional terms, their inclusion in the equa- 
tions does not influence the values of one-dimensional 
ones and the value of E Q . For the one-dimensional parts 
of the functions Sj, all terms in Eqs. 0DJ and (|H]1 
are of the same order, as it is easy to verify. Therefore, 
all they should be taken into account. But, in this case, 
(|3"5)) - (|4"l"j) is a complicated system of nonlinear integral 
equations difficult to be solved. Below, we will use the 
simplest zero approximations for s[ and Sg ■ 



s[ 1] (q x ) 



iV Nk lx Ci x (q x )2a 2 {-q x ) 
q x -2q x a 2 (q x ) 



(48) 



4 (q^qi) 



Relations (134 



W Nk lx Ci x (q x ) 2q lx a 2 (q x ) 

l-2a 2 (qj e (qi) + eo(q x +q x ) 

(49) 

and (US]) yield sf^qj ~ 



gi+9x v ?i qi+q-. 
(i) 



at g x < k rot {k rot is the roton momentum). So, S^ 1 takes 
values from ~ A^ 1 / 2 to ~ A 1 / 6 . 

It is easy to verify that the surface corrections A\ and 
A 2 have no effect on the solutions for Eq and a 2 (k): the 
values of A\ and ^2 by A 1 / 3 times less than those, which 
would render some effect. We expect that Aj>3 do not 
affect the solutions for aj>z- However, the boundaries 
change the "bulk part" of the equations for E and ay. 
the network of vectors k becomes by a factor of 2 d denser, 
and the Fourier-transform of the potential is multiplied 
by l/2 d . For this reason, the boundaries strongly influ- 
ence the values of E and dj. 



Besides WF (HD, WFs of the form (HI with * 



(and ki — > ki) are also the solutions. But such WFs 
are equal to zero inside the vessel, whereas ^0 (123) is 
nonzero. The theorem on nodes implies that the state 
with a WF without nodes corresponds to the minimal en- 
ergy. Therefore, the quantity ^0 is the exact and unique 
WF for the ground state of a Bose liquid and Bose gas 
in a box. 



III. GROUND STATE OF HE II: SOLUTIONS 



To verify the solution (|26l) , we will find the behavior of 
^0 in the case where the coordinate x of one of the atoms, 
for example the first one, varies from to L x . The sines 
standing before the exponential function in (|26[) set an in- 
homogeneous arc-like distribution of atoms in the vessel. 
But the density of He II is constant in the whole volume 
of the vessel except for a narrow strip (~ A) near the 
wall. The behavior of V^o should be analogous. Since the 



we need to consider the behavior of the ansatz ^ge^™' 

in (j2"6"1) . The function Sw (1^1, . . . , rjy) (l3"5"j) contains the 
one-particle dependence in each sum, and the values are 
of the same order. But taking account of the sums with 
<Sj>2 is a complicated task, and we neglect it. Then 

*o(a:i) ~ /(ri) = sin {k lx x x + 5 lx )e s ^\ (50) 



( 2 -) \ 
9x#0 VJV 



A 



+ 



+ E 

(2tt) 

+ E 

q x ,q y ,q^¥=o 



S[ 1} (q x +q z ) 



A 



(51) 



S- 1 (q,r + q„ + q,) ■ 



.) As was noted above, we have S[ 1] (q x + q + q 



S^in,. +%)/N 1 / 3 ~ S} l; (q x )/iV 2 / 3 . Therefore, three 
last sums in (jSTj) are of the order of the first one, but 
they depend additionally on y\ and z\. Let us consider 
only the first sum: 



(27T) 

si(r x ) =s x {x\) = 

q„ /L x 



(52) 



where = ±1, ±2, ±3, . . .. For (q x ), we use the zero 
approximation (|48|) . In view of the property (— q x ) = 
5^ (q^,) and since S^i^q^.) is real, the imaginary part of 
([52]) is cancelled, and we can change e l9x2:i — > cos (q x Xi). 
For investigating the qualitative behavior of the solu- 
tions, we choose the interatomic potential of two He 4 
atoms in the simple form 



U 3 (r) 



U b >0 r <a 
U bd <0 a<r<b 
r > b. 



(53) 



For the modern potentials [H, l34j . a « 2.6 A, Ubd 

II k . but the barrier and the well are more flat than 
the steps (|5"3"|) . Therefore, we will use in (|53p the lower 
values of a and \Ubd\- a — 2 A, t/jj = — 9K, and b = 4 A. 
The Fourier-transform of potential (1531 is 



47T 



(54) 



^3 (fc) = [C^fc — Ubd] [sin (afc) — a/c cos (afc)] + 

+ C/ M [sin (bk) - bk cos (bk)] . (55) 
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In the two-particle approximation, we have from ([43 



1 1 nv 3 (k)m 




(56) 



where ai = 4nmR 2 /2 d h 2 « 1.66/iO; B , d = 3, and i? « 
3.58 A is the mean interatomic distance for He II. For 02, 
we chose the sign "minus" before the root. At the sign 
"plus," we have nonphysical Eq = —00. 




Fig.l. Values of the functions F = f(x) (o) and F = s\{x) (★) by 
formulas l(50)l . (|52T> for potentials (14) 1 and ((53j) at a = 2 A, Ub = 
1600 K, U b d = OK, U s — > 00 and the vessel size L x = 10 4 R. Values 
of f(x) and si(x) at x/L x £ [0.5,1] are obtained by the reflection 
of the curves in an upright mirror, positioned at x = 0.5L X . 

The results for the functions f{xx) and Si(xi) (I5TJ1) . 
(l52"j) at L x — 1{) A R and several potentials are given in 
Figs. 1 and 2. The values of the functions at L x = 10 5 R 
and L x — 10 3 R differ from the drawn ones very slightly 
(by < 0.1%). So, the curves in the figures are true for 
all L x > 10 3 .R. As seen from Fig. 1, f(xi) is equal 
to 0.5 everywhere, except for a narrow region near the 
wall. While approaching the wall, f(x\) oscillates with 
constant period and increasing amplitude and becomes 
zero on the wall. The dependence on the potential is as 
follows. An increase of the interatomic barrier Ub causes 
the growth of the amplitude of oscillations /(xi). The 
well depth Um does not affect the results, if \Um\ *C Ub- 
As the distance to the wall increases, /(xi) approaches 
the constant: /(xi) — 0.5. This holds with a high accu- 
racy increasing with L x . For example, at L x = 10 4 i?, 
the deviation from 0.5 occurs only in the sixth decimal 
point. As x\ —¥ L x , the function /(xi) behaves itself as 
at x\ — > 0. 




Fig. 2. Values of the function f(x) J50J, J52J near the wall (x = 0) 
at L x = 10 4 R, Us — > 00 and potential l53t with a = 2 A, b = 4 A, 
U bd = -9K, and three barriers: U b = 800 K (o), U b = 2000 K (*), 
and U b = 2900 K (o). At x > 6 A, curves 000 and * * * almost 
merge. Values of x are given in A. 

It is of interest that f(x\) oscillates near the wall. This 
is because the short-range order, and the period is pro- 
portional to a/R. Large oscillations at Ub = 2900 K are, 
most likely, nonphysical and are related to the neglect 
of corrections, which are large at such a potential. This 
reflects the behavior of the density p of helium, since 
p{x\) ~ j dyidzidr2 ■ ■ ■ dr o\ 2 = const x f 2 (xi) (the 
behavior of p{x\) near the wall is also affected by the 

integral of e Sb , but we did not study it). Similar oscilla- 
tions were obtained in [2l| , while modeling the properties 
of a film of helium adsorbed on the substrate. 

We obtained the solutions also at finite U s . At 
U s > 1000 K, the solutions are practically invariable with 
increase in U s , except for the values on a wall, where / 7^ 
0. As U s increases, the function s\{x) approaches a con- 
stant (~ 10) at the point x — 0, and sin^izXi + 8\ x ) — > 
0. As a result, we have f(x — > 0)\u s ^oo — > 0. With re- 
gard for the higher corrections, such a behavior will be 
conserved: as U s — > 00, we have <5 — ?► 0, such S affect 
slightly C\{q) and all results, because 5 enter the equa- 
tions through the function C\ (q) . In addition, the terms 
with large q are cut in sums by a,j. Due to such proper- 
ties, the transition to opaque walls (U s — > 00, 6 — > 0) is 
correct and makes ^0 to be arbitrarily close to zero on 
the walls. 

Thus, though the bare function sin (k\ x xi + 6i x ) is far 

from to be a constant, the account for S{ (q x ) makes the 
theoretical density of atoms to be constant everywhere, 
except for a narrow (~ A) region near walls. Therefore, 
sin (kixXi + Si x ) is a "good" bare function, despite the 
anxiety (9rlTl|. It satisfies the equation for free particles. 
If the interatomic interaction is "switched-off," ^0 is re- 
duced to product (|20]) of such functions, as it should be. 
Apparently, this is the best bare function: at the other 
bare functions, we obtain the same solution for ^q, but 
in other form. However, the properties of the solution are 
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not improved in this case, whereas the equations for Sj 1 ^ 
become too complicated. The circumstance that the cor- 
rection s[ 1 \q x ) "improve" the behavior of f(x\) seems 
natural, because the Schrodinger equation (|T8]) involves 
the interaction of atoms and must make their distribution 
uniform, but s[ (Qx) is generated just by the Schrodinger 
equation. 

We can find E (|38]>. (|42|) with regard for A x « and 
a,2 ([56"]) . However, the potential of a semitransparent ball 
([53"]) is "bad" for this purpose in view of the asymptotics 
vs(k — > oo) ~ cos (ak)/k 2 , at which the value of E$ (|4"2l. 
([56]) is uncertain. For 3D-systems, the suitable asimp- 
totics is v s (k —> oo) ~ l/k p ,p > 3. Therefore, we will 
not estimate Eq here. 



IV. EXCITED STATE OF HE II 

We now find the WF of helium II in a vessel for the 
state with one phonon. Under periodic BCs, the exact 
solution for the WF takes the form [23 



* k (ri,...,r^) = k *c 



(57) 



(2tt) 



6 2 (k 2 ,k) 



k 2 #0,-k 



N 



pk 2 p-k 2 -k+ (58) 



(2tt) k 2 +k 3 +k#0 

E 

k 2 ,k 3 #0 



fr 3 (k 2 ,k 3 ,k) 



Pk 2 Pk 3 P-k 2 -k 3 -k + 



The approximate 0k of such a type, which consists of the 
zero approximation 

(59) 



P-k 



and a simplified correction with 62, was first proposed by 
R. Feynman [U, H| . The arguments to favor of solution 
([ST)) . ([55]) are as follows: 1) it satisfies the Schrodinger 
equation for N interacting Bose particles and 2) it is an 

eigenfunction of the momentum operator P = —ih^ V rj 

j 

with the eigenvalue Hk; 3) at the "switching-off" of the 
interaction, it is reduced to WF ([5U1) describing N free 
Bose particles, from which N — I particles are immovable, 
and a single one moves with the momentum k; 4) an 
analogous solution is obtained in the operator approach 
3; 5) the completeness [l7| of collections (TK)1) implies 
that the solution for 0k should have only such a form. 
However, the structure of (|57[) is only the assumption, 
though no solution with another structure was found. 
The solution for -0k with "shadow" variables [IH, is 
equivalent to ([57]). (1551) . 

For the fluid in a vessel, the WF "Jk should be sought 
anew, since *k (JSjJ, (|58[) is not the solution of the 
Schrodinger equation. 

The simplest way to find solution (|57l) . (|58p is to as- 
sume the structure of (|57[) and to substitute the WF of 



free particles ((59)) . as the bare 0k, in the Schrodinger 
equation for interacting particles. The equation prompts 
the form of corrections to 0k- 

The reasoning for particles in a vessel is analogous. 
We start from ([57]) with f ([26]). The study of different 
possibilities has shown that "0k must be sought in the 
form 

8 N 

-0k = —j= ^ cos (K x Xj) cos {K y y.j) cos {k rz z ) + <50 k 

V 3=1 

(60) 

with the unknown <50k- In this case, k r = (fc rx , k r , k r ^) 
takes values from series ([3]), since we must obtain the 
solution for free particles in a box with law ([3]) at the 
switching-off of the interaction. We note that 

8 N 

cos (k ric Xj ) cos (k Ty yj ) cos (k Tz Zj ) = 



(61) 



j o_k r + 7 permutations, 



where the permutation means p_k with the different sign 
of one or several components of the vector k r . We now 
represent 0k (|60[) in the form 

(62) 



0k = 0k + 7 permutations, 



0£ - P-k + 6i>l (63) 

It is easy to see that solution (|62p , (1531 describes a three- 
dimensional standing wave decaying into eight counter 
running waves. 

Substituting WF ([37]) in the Schrodinger equation with 
Hamiltonian (]12p . we obtain the equation for 0k: 



^- {A,0 k + 2(V,0 k ) x VjS} = i?(k)0 k , (64) 



S = + S b + S , 



50 = 1x1*/ 



(65) 
(66) 



where -E(k) = E — NEq is the energy of a quasiparticle. 
-B(k) is obviously independent of the signs of components 
of k. Therefore, Eq. ([M]) is separated into 8 equations: 
for 0£ and 7 permutations with the same energy E(k). 

The analysis indicates that the solution for 0^ has the 
following form: 

k =0£ + &o(k)+ Qi(q,k)P-q-k + 

q#0,-k 

q+ qi +k#0 

<?2(q,qi,k) 



E 

q,qi#o 
q+qi+q 2 +k#o 

E 

q,qi,q 2 #o 



iV 



~PqiP— q-qi-k 



(67) 



Q3(q,qi,q 2 ,k) 

N 



"Pit P12 P-i-ii -q 2 -k 
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where ip^ is the solution (|58[) for periodic BCs (but with 
kj running series ©), &o and Qi are the corrections from 
boundaries. Here, q is quantized according to (|2j (like 
2nj/L), and q - and k are quantized by ^ (like wj/L). 



With regard for the properties of p k , we can directly 
verify that <f k ([ST)). ([62]). (lrJ7|) is the exact solution of the 
Schrodinger equation, if the functions E(k), bj, and Qi 
satisfy the equations 



6o(k)e(fc) --2fc 2 ^ 1) (-k)<5 k , k e 



2ki x k x iCi x (-k x )V NS k ^ - 
^g 2 Q 2 (-k,q,k)^ k e 



iV 

(2tt) 



q#0 



- J2 2(k + q) 2 g 1 (q,k)5i 1) (-k-q)^ k , k . + 

q^O ,-k 
(2*0 _ 

+ ^2 2VNk lx q x iC lx (q x )Q 1 (-k - q x , k)c$ k , k .= + 
+ (x->y,z), (68) 



(*) 



e(fc) = e (fc) ~ ^ E &2(k 2 , k)2k 2 (k 2 + k) - 

k 2 #0 

, W 

- -$>fc 2 & 3 (k 2 ,-k 2 ,k)- 

(2*0 

- E 2k lx (k x + q x )iC lx {-q x )Q\(q_ x , k) - 

2 (2T) 

- 7 =^Q 1 (q,k)(k + q)x 

x [qSf 5 (-q) + 2(k + q)Sf >(-q, ~k)] + 

(2*0 

+ E 4fci x ? :r iC'i x (g' x )Q2(-q :E) q :c ,k) - 



(2*0 



(69) 



7 =E^(-q,q,k)^ 1) (q) + ( 



Qi(q,k) [e(fc)-e (k + q)] = 
2iki x C\ x (q x )5< l ^ [-k x + 2g a; 6 2 (q,k)] 

2S< 1) (q) 



4fc 2 



[qk-2g 2 6 2 (q,k)]-^=4 1) (q,-q-k) 
W v jV 



(*0 



[Q 2 (q, qi ,k) + Q 1 (q + q 1 ,k)ViV5 1 (1) (-q 1 )] + 
6g 2 Q3(q, q l7 -q!,k) + 

4<Z 2 Q 2 (q - qj, qi.kjVjVSPfqj) + Qi(q + q x ,k) x 
4( Ql + q + kfVFsf ^-q^ -q - k)} 

(2*) 

2fci a! «<7i !1 .(p a .) [(-A; x - g x +Px)Qi(q.~ P^jk) + 
2p a Q2(q- Px,P x! k)] + (a;-> j/,z), (70) 

6 2 (k 2 , k) [e(fc) - e (k 2 ) - e (k + k 2 )] = 
kk 2 a 2 (k 2 ) - k(k + k 2 )a 2 (k + k 2 ) - fc 2 a 3 (k, k 2 ) - 



— £ {2 qi ( qi +q + k) 



qi#o 



(*) 



3k 3 (k 2 +k 3 +k) 



iV 



6 3 (k 2 ,k 3 ,k) 



+ Qi(k 3 ,k)(k 3 + k) ^(ka-kaJS^C-ka.ka) + 

+ 3(k 3 + k)^ 1) (-k 3 ,k 2 ,-k 2 -k)] + 

+ 2Q 2 (-k 3 ,k 2 +k 3 ,k)(k 2 +k 3 ) x 

x k 3 ^ 1) (k 3 )+2(k 2 +k 3 )^ 1) (k 3 ,k 2 )] + 

+ 3fc 2 sf ) (k 3 )Q 3 (-k 3 ,k 2 ,-k 2 -k,k)} + 

(2*) 

+ ^2 2k lx iC lx (p x ) [3p x Q 3 (-p a ,,k 2 ,-k 2 -k,k) + 

+ 2(fc 2 ,+ P;c )g 2 (-p :c ,k 2 + p ;c ,k)] + 

+ higher corrections + (x — > y, z), (71) 

where e(k) = 2mE(k)/h 2 , k e is the wave vector with 
even components (i.e., they are multiple to 2w/L), q x = 
q x i x ,q y = q y i y ,q z = q 7 \z, and analogously for p x , p y , 
and p z . By (x — > y,z), we denote the same terms as that 
with the separated ^-component, but with the changes 
x — > y and x — > z. All terms of Eqs. (|6"8]) - ([7"T1) . except 
for the first one, contain a product of two wave vectors 
(for example, q 2 , k lx k x , or (k 3 + k)(k 3 — k 2 )). These 
wave vectors must be nonzero. 

Equations (I68l) - (|71[) are written in the approximation 
of "two sums in the wave vector" , at which the series con- 
tain the functions a 2 , a 3 , fe 2 , & 3 , Sj]) 3 , and Qz< 3 and do 



not include the corrections Oj>4, fej>4, <5j> 4 , and Qi>4- 
Taking these corrections into account will change the 
equation for 6 2 , by adding the correcting sums with &4, 

$4 , and Q4. 
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If the interaction is switched-off, ^k reveals a compli- 
cated structure, which should be reduced to the solution 
for free particles; but we omit this point. 

Let us consider the properties of Eqs. (|57)) -([7T |) . Solu- 
tion (|58p was obtained under periodic BCs and contains 
only the functions 6j>2- Solution (|67[) corresponds to the 
zero BCs and contains the additional functions 60 and 
Qi. In view of the "one-dimensional" form (|46p . {4"?]) 
of the functions S^, relations ([gg j) -([TT ]l yield that the 
functions Qi can be also considered "one-dimensional." 
In particular, 



Qi(q,k) ps Q!(q,k)((5 q:qx 



<W + <W)- 



(72) 



The "two-dimensional" Q\ with q of the form (q x ,q y ,0) 
and the "three-dimensional" ones with q = (q x ,q y ,qz) 
are less than the one-dimensional ones by ~ N 1 / 3 and ~ 
N 2 / 3 times, respectively. Due to the smallness of nonone- 
dimensional Q%, their inclusion in the equations does not 
affect the values of one-dimensional ones. The estimates 
indicate that Qi is less than k\k and affects slightly e(fc), 
since the bulk corrections (the second and third terms on 
the right-hand side of (|6U]) ]) are of the order of eo(k) ~ 
kk rot , i.e., by ~ N 1 / 3 times larger. In Eq. (jTTj) , the 
surface corrections for & 2 are also less than the bulk ones 
by ~ N 1 / 3 times. The smallness of surface corrections Ai 
and Qi is obviously understandable and is related to the 
thinness of a near-surface fluid layer (Fig. 2) as compared 
with the sizes of a system. 

Therefore, we can set Qi = in Eqs. ([57 |l -([7I ]l . Then 
we obtain the well-known Vakarchuk-Yukhnovskii equa- 
tions [27] for e(k) and bj. However, the wave vectors kj 
in these equations should be quantized by the law tt/L 
Q instead of the usually used 27r/L-law ©, which is 
valid for the periodic BCs. In the zero approximation for 
^0 and ^k, we set dj>3 = and bj = in the equations. 
Then we obtain the Bogolyubov formula for the energy 
of quasiparticles, 



E(k) 



fh 2 k 2 



\ 2m 



nv 3 {k) h 2 k 2 
2 d m 



(73) 



where, however, the potential has the additional factor 
l/2 d due to effect of the boundaries (in the general case, d 
is the number of nonperiodic coordinates). This formula 
is true at a weak interaction. Comparing formula (|73[) 
at potential (|53|) with the dispersion curve of He II as 
k — >• 0, we obtain Ub — 283 K, whereas we have Ub — 
35 K without the factor l/2 d . The first estimate is more 
plausible. 

The Bogolyubov-Zubarev field model Q gives the same 
results: Since the potential should be expanded in Q in 
the Fourier series (TTBI . we replace v 3 (k) — > v 3 (k)/2 d and 

(2*0 (tt) 

—> m formulas |3j and obtain the formulas for 
E(k) {73]) and E @2J with 03 (ft) {SB]). 

For a Bose gas in a trap, the sound velocity was mea- 
sured 35]. Its value agrees approximately with the for- 
mula c ~ yj nv 3 {Q) /2m. We cannot compare these results 



with ours, because our calculations did not involve exter- 
nal fields, whereas the field of a trap plays a significant 
role in the experiment in (35[. 

Relations d(H)/dv 3 (k) JdH/dis 3 (k)), 
S(k) = (Pkp-k) and (HID yield 27] S(k) = 
1 + (2 d+1 /n){dE /dv 3 (k)). Using this result, {38]), 
and (|42l) . we will find the connection between az (k) 
and the structural factor S(k). In the zero approx- 
imation, it is the same like that for periodic BCs: 
1 — 2a2(k) ps 1/S(k). With regard for this relation and 
(1561) . formula (|73[) is transformed into the Feynman 
formula 



E(k) ps H 2 k 2 /2mS(k), 



(74) 



which agrees approximately with the experiment. 

If one particle is in the box, k is restricted from above 
by the value k max = yfd>mU 8 /h (see {25]) ]). For He 4 atoms 
at U s = 25 K, we have k max ps 3.6 A" 1 . So, the observed 
break of the dispersion curve of He II at k ps 3.6 A can 
be caused, in principle, by that a phonon must have also 
some k max at a finite barrier of the wall, rather than by 
the decay mechanism [36| . 

Additional comments for ifi^ can be found in Appendix 

A. 



V. ROLE OF BOUNDARIES 

Let us return to the effect of boundaries, which is deli- 
cate. As usual, the works on quantum fluids and crystals 
do not consider the boundaries and use periodic BCs. 
The discussion of a role of boundaries can be found only 
in a few works. In particular, it was noted by J. Ziman 
[37] ] that there exists a bulky proof of the correctness of 
periodic BCs (but it is not presented), and that the iden- 
tical density of states (for quasiparticles) under real and 
periodic BCs testifies to the applicability of the latter. 
K. Huang [11] gave the approximate reasoning in favor 
of periodic BCs, and expressed the cautious hope for that 
the periodic BCs do not lead to erroneous results. 

It was mathematically shown above that this hope is 
not justified: the boundaries change significantly the so- 
lutions. How we can imagine this? According to Sec. 2, 
BCs reveal themselves at the expansion of the potential 
EMI r i ~ r 2|) m a Fourier series: the expansion is real- 
ized with wave vectors multiple to 2tt/L under periodic 
BCs and to tt/L under zero BCs. This becomes clear in 
the following example. Let us consider oscillations of two 
identical strings with length L: one string is unclosed and 
has fixed ends, whereas another one is closed. The un- 
closed string admits the standing waves with the integer 
number of the waves (L = = 1,2,3...) or half- 

waves (L = j\/2,j = 1, 3, 5 . . .) on the string. But, for 
the closed string, only waves of the first type are realized. 
Periodic BCs and zero BCs correspond to closed and un- 
closed strings, respectively. We see that the collections of 
natural oscillatory harmonics are different in both cases. 
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The ground state of a Bose fluid can be compared with 
an oscillating string, since "Jo — ex P (2 a 2(k)/OkP-k/2) 

k 

has the structure of the WF of the ground state of oscil- 
lators with "coordinates" pk- The boundaries change the 
collection of harmonics of the GS. Due to the interaction 
of harmonics, their intensity is changed as well. There- 
fore, the solutions for \l/o are different for periodic and 
zero BCs. Moreover, due to the coupling of the equations 
for the ground and excited states, the harmonics of the 
GS affect the excited states. Thus, though the density 
of excited states is the same for periodic and zero BCs, 
the sets of k do not coincide and the dispersion curves 
are different. As is seen from the comparison with closed 
and unclosed strings, we deal with a topological effect. 

It is usually considered that the boundaries must not 
affect the bulk properties of a system, since the volume of 
a near-boundary region is negligible relative to the total 
volume of a system. The smallness of the near-boundary 
region yields, in particular, the dimensional formulas of 
the form (3.lV-(3.3) from [38| . in which the bulk charac- 
teristic (in 38], it is the ground-state energy per atom) 
of the system is invariable as N, V — > oo, N/V — const. 
This implies that the shape of boundaries does not affect 
the bulk properties. We agree with this, it corresponds 
to the smallness of surface corrections Aj and Qj in the 
above-obtained equations for ^ and Vf^. Since any infi- 
nite system (N, V = oo) has no boundaries, we make the 
second conclusion: the very presence of boundaries does 
not affect the bulk properties. One would think that this 
conclusion follows from the proofs [3{| of the existence 
of a thermodynamic limit too. But we meet here the ab- 
sence of a full rigor. A system with boundaries differs 
from an infinite system topologically: the latter has no 
boundaries and is topologically closer to a closed system. 
If the transition N, V — > oo means that to an infinite sys- 
tem, it is implicitly assumed that a change in the topol- 
ogy does not change the class of functions describing the 
bulk properties of a system. This should be proper for 
quantum systems at high T: the physical processes for 
them are conditioned by atoms, and an atom after the 
collision with a wall collides with other atoms and loses 
the memory of the wall. However, as T — > 0, the physics 
of systems is usually conditioned by waves (collective os- 
cillations). The standing waves are modulated by walls 
and, therefore, "remember" the walls. Because the col- 
lections of harmonics for waves in an open and closed 
system are different, the bulk properties of these systems 
can be different as well. The above-presented calculation 
for He II confirms this. Since we deal with a topological 
effect, it should be valid not only for quantum fluids, but 
also for gases and crystals. In our opinion, the proofs of 
the existence of a thermodynamic limit mean only that it 
is possible to pass to the limit N, V — > oo, N/V — const 
in the formulas. But, in this case, we remain in the class 
of functions obtained for a finite system. If the infinite 
system is described by other functions, such a transition 
does not allow one to describe the infinite system, i.e., 
this transition is incorrect. Of course, such a problem 



does not arise for cyclic systems, because the transition 
A, V — > oo does not change their topology. The main 
point consists in that the effect is related to the topology 
of a system as a whole, rather than to the influence of 
a thin layer of near-surface atoms. In this meaning, the 
effect is a bulk one. 

The role of BCs is also discussed in Appendix. 



VI. THERMODYNAMICS 

Under periodic BCs, the thermodynamic quantities are 
calculated from the free energy [I(J 

F = -k B Tj\n(l + n q ) — . (75) 

Under zero BCs, law ^ holds. Therefore, the quan- 
tity 2ir in ff75|) should be replaced by ir, which increases 
the value of the integral by eight times. However, if we 
change the sign of one or several components of k in the 
WF of a standing phonon, the WF is not changed. In 
other words, eight states with different signs of the com- 
ponents of k are equivalent. Therefore, we must integrate 
only over the sector k x ,k y ,k z > 0, i.e., integral (1751 de- 
creases by eight times. We arrive at formula (|75|) and the 
known formulas for entropy, heat capacity, etc. 



VII. POSSIBLE EXPERIMENTAL TESTS 

Formula (|73[) implies that the dispersion curve of He 
II must be different for periodic and nonperiodic bound- 
aries. In particular, it should be different for a mono- 
atomic film of He II on a plane surface (d = 2 in Eq. 
(|73|l). on the lateral surface of a cylinder (d — 1), and 
on the surface of a torus (d — 0). Therefore, the heat 
capacity of these films at the same temperature must be 
significantly different. A more interesting effect is also 
possible. It is necessary to form a monoatomic film of 
He II on the surface of a torus and then to cut the torus 
(similarly to the breaking of a ring). As a result, some 
boundaries must arise. In this case, the system passes 
from the state with d = to that with d = 1. There- 
fore, the energies of phonons and of the ground state are 
sharply changed, the ensemble of phonons will be rear- 
ranged, and the temperature of the system will jump by 
a value comparable with the initial T. If we restore the 
cut torus, the inverse processes will run. This is a phase 
transition that can be called a topological one. Its ob- 
servation is quite real and can be a test for the effect of 
boundaries. 

Our analysis shows that if the film thickness is more 
than several atomic layers, the properties of a helium 
film on the cylinder surface must be the same as those of 
helium in the cylinder bulk. 



12 



VIII. DISCUSSIONS 

Of interest is the following property. The ground state 
is described by WF (|26]) with the product of N triple 
sines with the wave vector ki, and the excited state is 
presented by the same WF with the factor ifa ■ The sines 
of the GS set a standing wave in the probability field. In 
the zero approximation, ip-^ is equal to p_k with permu- 
tations and "is convolved" with the sines of the GS, by 
resulting in a superposition of states, whose structure co- 
incides with that of with the difference that the wave 
vector for some sines differs from ki. We may roughly 
say that the excitation is simply the replacement of one 
of the standing waves of the GS by a wave with larger k. 
Such a structure of WFs means that the GS is formed 
by N identical standing "phonons" with the wave vector 
ki. In the excited states, a part of phonons is replaced 
by phonons with different k. Therefore, we may say that 
the ground state and excited ones of the Bose fluid in a 
rectangular vessel contain N harmonics (phonons), and 
the difference consists in the composition of harmonics. 
In this case, the meaning of the notion of "excitation" 
is changed, because the GS becomes, in a certain sense, 
"excited" , since it contains N harmonics with minimal 
k. 

With regard for corrections, the harmonics of the GS 
are not identical by structure to the harmonics of weakly 
excited states: the corrections for harmonics of the GS 
are collected in the exponent, whereas this is not the case 
for the excited states; and the equations for corrections 
are different. This is not amusing, because the GS is the 
resonant state with N identical interacting "phonons". 
Their structure must be different from the structure of 
single "ordinary" phonons described without regard for 
their interaction. 

By passing through helium, neutrons create phonons. 
This process occurs with the conservation of the momen- 
tum. According to our solution, phonons are standing 
waves without definite momentum. However, the WF of 
a phonon is a sum of eight traveling waves, each with 
a definite momentum. Neutrons create, apparently, just 
such running waves. In a small region of helium, we see 
a lot of quasiparticles in the form of moving localized 
wave packets with a certain momentum. But if we con- 
sider helium on the whole, these quasiparticles turn out 
to be parts of the total system of standing waves filling 
the vessel. 

We now return to the method. As was shown (l7l |27j. 
the collections p\ !Lj (|19p are functionally independent, if 
kj are quantized by the 27r/L-law @. But, in our sums, 
the wave vectors are quantized sometimes by the 7r/L-law 
(O; in this case, the collections p^, are not func- 
tionally independent. For example, pk with odd k (with 
components of the form (2j + l)w/L) can be expanded in 
collections (19) with even k (with components of the form 
2jw/L). But, in this case, we are faced with difficulties 
indicated below. Therefore, we find a solution without 
expansion. Such a means is quite correct. Indeed, let us 



seek a function in the form 

F(x) = Y / A 3 f J (x), (76) 

3 

and let the Schrodinger equation is reduced for it to the 
form 

Y / B l ({A J })f l (x) = 0, (77) 
i 

and we determine coefficients Aj such that Bi = for 
all I. Then it is obvious that F(x) is a solution irrespec- 
tive of whether the functions fj{x) are independent and 
whether they form a full collection. We have acted just 
in this manner. As usual, the full collection is used in 
order that all solutions be found. But, under a certain 
structure of the equation, it is easier to find solutions if 
fj(x) have some particular structure and are dependent. 
It seems that this is our case. We have used the expan- 
sion in the variables p^ j p^|) with k ~ ttj'/Z (0). Part 
of them, namely pk (|19p with even k, forms the full col- 
lection, and the other variables are additional. Since the 
collection of functions, in which we expand, contains the 
full collection, no solutions should be missed. 

We can try to use the independent basis functions as 
follows. Let us expand the potential in a series 

I E U 3 (\ri -Tj\) = a<pi + ^c 2 (k)<p 2 + 

ki+k 2 #0 

+ C3(kl,k 2 )¥>3 + ... (78) 

ki,k 2 #0 

with k, ([2]), where the functions [13] 

tpi=Ax, ip 2 = A 2 (pkP-k - 1), (79) 

^ = tI (/N pk i pk * p -*- i « + 2 ~ (80) 

— (PkiP-ki + Pk 2 P-k 2 + Pki+k 2 P-ki-k 2 )) i • • • 

[Aj = const) form a complete orthonormalized basis for 
the expansion of any translation-invariant Bose symmet- 
ric function of the variables r 1; . . . , rjy. However, such a 
method turns out not the best one: due to the presence 
of divergent terms, and because many basis functions in- 
clude one- type terms of the form pkP-k, e ^c. such that 
it is difficult to gather and to evaluate them. In ad- 
dition, the distinction of the zero and periodic BCs is 
insignificant in such an approach: for example, the co- 
efficient C2 differs only by a small relative value ~ o>/L, 
where a is the interaction radius. For He II, we have 
a/L ~ l/N. Therefore, in order to reveal the distinction 
of the BCs, we should consider ~ N corrections in ex- 
pansion (I78p . which is impossible to perform. In other 
words, the language of the functions ipj is not suitable 
for taking boundaries into account. 
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Note also the following. The expansion in the 2ir/L- 
harmonics reproduces exactly the potential at the zero 
BCs at \xi - Xj\, \yi - yj\, \zi - z 3 \ < L/2. To specify 
the potential in the domain > L/2, we must additionally 
attract the odd 7r/L-harmonics from (fT5|) . In this case, 
the function is set by different collections of harmonics 
on different intervals. Therefore, the description with the 
use of only even 27r/L-harmonics in the whole domain is 
not efficient and leads to the difficulties related to series 
(|78| . As opposite, the expansion ((T5)) reproduces the 
potential in the entire domain on the basis of the single 
collection of harmonics and, therefore, is better. 

IX. CONCLUSIONS 

We have studied the influence of the walls of a vessel 
on the microstructure of He II and obtained a new class 
of exact solutions for the wave functions describing N 
interacting Bose particles in a box. 

It is shown that the explicit account for boundaries 
changes the equations for the Bose liquid and the calcu- 
lated values of the observed bulk quantities. Thus, the 
models with periodic boundary conditions should be re- 
considered. 

Probably, the wave function of the ground state of Bose 
crystal with rectangular lattice can be constructed [4l| 
on the basis of the wave function ([26| of Bose liquid. 

Recently, we have solved the problem in an approach 
with independent basic functions [421 ] . In the limit of 
a weak interaction, the solution coincides with that ob- 
tained above. The further studies will help us to better 
understand how the boundaries affect the properties of 
Bose systems. 

The author thanks N. Iorgov, K. Szalewicz and S. 
Yushchenko for valuable discussions. 



X. APPENDIX A 

It follows from (|67|) that the functions Q;> 2 satisfy the 
symmetry relations 

Q 2 (q,qi,k) = Q 2 (q, -qi - q- k,k), 
<23(q,qi,q 2J k) = <23(q,q 2 ,qi>k) = 

= Q3(q, -qi - q 2 - q - k, q 2 , k) = 
= Q3(q,qi,-qi - q 2 - q- k,k), 

and the analogous relations can be written for higher Qi. 

Let us consider b (k) (|68|) . As seen from (|68"1) . b (k) 
is nonzero, if all components of the vector k are multi- 
ple to 2tt/L. It follows from the equations that bo does 
not affect the dispersion curve e(k). Since e(k) ~ kk ro t, 
relation (|68]) implies that, for the one-dimensional k, we 
have 60 ~ N 1 / 6 at k ~ k rot . For the two- and three- 
dimensional k, bo is less, respectively, by N 1 / 3 and TV 2 / 3 



times. For the irregular configurations of atoms, typical 
for a fluid, p^ ~ 1. Therefore, for the one-dimensional k, 
we have 60 (k) >• pk, i-e., bo determines the quantity ipk 
([rJT|). At the normalization, bo "hides" itself, and we ob- 
tain ~ 1- By the structure of ipk, such states are close 
to the ground one, but their energies are different at dif- 
ferent k and larger than Eq. The nature of this anomaly 
is not clear. It can turn out fictitious, since Vf^ (1571) . (pO)l 
with different k are not completely orthogonal due to 
the corrections with Q[. These corrections cause a small 
mixing of different states: / dil^k^^^ ~ N~ 1/3 at 

In addition to the regular solution, Eqs. (|M )) - ([7T|) have 
a huge number of pole solutions related to the nullifica- 
tion of the denominators of the functions 6 2 , 63, and Qi 
in (|69|) . For example, let we have a set of pole solutions 
of the form e(fc) = e (k 2 ) + e (k + k 2 ) +5e(k, k 2 )/N for a 
given k with different k 2 and with Se ~ 1. It can be veri- 
fied that these solutions disappear in the thermodynamic 
limit, but they form a new solution class for a system 
with finite size. However, these solutions are nonphysi- 
cal. Indeed, we have 6 2 (k 2 ,k) ~ N for pole k 2 . Hence, 
the term 62( ^ k) Pk 2 P-k 2 -k ~ \/iVpk 2 P-k 2 -k, which is 
greater than the term p_k by ~ y/N times, becomes the 
main term (~ \/~N) in expansion (|58j) . The quantity i/'k 
corresponds to distorted two-phonon state. The same is 
true for poles of Q 2 . The poles of 63 correspond to dis- 
torted three-phonon states. Moreover, the poles of Qi 
correspond to a phonon with a different wave vector. 

XL APPENDIX B. MODEL OF POINT 
PARTICLES 

The explicit account for boundaries was performed 
within the exactly solvable one-dimensional model (l2l . 
[l3| for bosons with the point interaction (point particles 
with the interatomic potential as the (^-function). For 
such a system, the periodic BCs turn out, in essence, 
to be equivalent to the zero BCs: the values of Eo for 
both BCs differ from each other only by a negligible value 
~ Eq/N. At the same time, it was obtained above for 
the Bose fluid that the values of Eq differ strongly under 
periodic and zero BCs. 

Such a difference may be caused by the difference of in- 
teratomic potentials. In our opinion, the principal point 
consists in the following. Consider a cycli c sy stem. The 
solution for the WFs of point particles is [38] : 

N 

^ — » i y~* kp. xi 

* (fe) (x) = ^a(P)Pe .=1 , (81) 
p 

where P are all permutations of k\. WF (|5Tj) is a super- 
position of states of N free particles. The energy, 

E={h 2 /2m)(k 2 1 +kl + ... + k 2 N ), (82) 

is the same for all these states. At the nonpoint inter- 
action, the WF of the ground state is given by (@|, (JSJ. 
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If we expand this WF in a Fourier series, we obtain a 
superposition of states (in ID) 

JV 

* o (a0 = ^a(P({x l })) ]T A(k u ...,k N )e^ klX \ 

P ki...kjsr 

(83) 

with ki = 2irl/L and energy which is different for 

different terms. This "J/o cannot be presented in the form 
of a superposition of plane waves with the same energy. 

In the presence of walls, a WF of point particles was 
constructed [l2| by the optical analogy as a sum of plane 
waves (f8Tj) and their all reflections with the same en- 
ergy. Therefore, a wall does not change the energy of 
the system. According to [l2j|, the solutions for N par- 
ticles in a vessel with length L coincide with definite so- 
lutions for 2N particles in a periodic vessel 2L in length. 
However, the waves in (|53")) have different energies, and 
the addition of reflected waves will lead to the redistri- 
bution of terms and to a change in the energy of the 
system. This can be seen at another angle of vision, 
if we take into account that, according to ((H) and ([S]), 
V&o = exp (J2 a2(k) j OkP-k/2) is the WF of the ground 

k 

state of oscillators with coordinates p^, and the correc- 
tions Oj>3 are related to the interaction of oscillators. Let 
us have a pair of oscillators with the same frequency uj\ . 
If we add to them a pair of the same oscillators, which 
does not interact with the first pair, then the frequency 
of all oscillators remains to be This corresponds to 



the solution for a point interaction. But if we include the 
interaction between pairs, then the eigenfrcquencies of 
the system are changed, as is known from mechanics. In- 
stead of two frequencies, we will have four ones, or more. 
This case corresponds to a nonpoint interaction. 

The absence of the effect of boundaries under a point 
interaction can be imagine, probably, in the following 
way: The effect is related to a certain game of harmonics; 
however, the point interaction has no dispersion {v{k) = 
const) , and no game is present. 

We have tried to qualitatively explain the effect. But 
these analogies are rather rough, and only the exact so- 
lution of the problem can be a proof. 

Since \l/o of the systems with point and nonpoint in- 
teractions are different, the solutions for Eo and E(k) 
should be different too. However, it turns out that the 
solutions [HI, EH, 113 for a weak point interaction are close 
to Bogolyubov ones. This is related to the fact that Bo- 
golyubov approach does not use WFs, and it is possible 
to describe the systems with point and nonpoint inter- 
actions in the same way (in the first case, the potential 
can be considered in the r-space not passing to Fourier 
transforms). Therefore, it is of importance to find the 
solutions for a noncyclic nonpoint system namely in Bo- 
golyubov approach and to determine to which solutions 
they correspond: to the traditional solutions for a cyclic 
system or to the above-obtained ones. We now study this 
problem. 



L. Landau, J. Phys. USSR 5, 71 (1941); 11, 91 (1947). 
N.N. Bogoliubov, J. Phys. USSR 11, 23 (1947). 
N.N. Bogoliubov and D.N. Zubarev, Sov. Phys. JETP 1, 
83 (1956). 

R. Feynman, Phys. Rev. 94, 262 (1954). 

R.P. Feynman and M. Cohen, Phys. Rev. 102, 1189 

(1956). 

H.R. Clyde and E.C. Svensson, in Methods of Experimen- 
tal Physics, vol. 23, part B, edited by D.L. Price and K. 
Skold (Academic Press, New York, 1987), p. 303. 
L. Reatto, J. Low Temp. Phys. 87, 375 (1992). 
MD. Tomchenko, arXiv:cond-mat/0904.4434. 
A. Bijl, Physica 7, 869 (1940). 
MD. Girardeau, J. Math. Phys. 6, 1083 (1965). 
K. Huang, Statistical Mechanics (Wiley, New York, 
1987), Chapter 19. 

M. Gaudin, Phys. Rev. A 4, 386 (1971). 

M. Gaudin, La Fonction d'Onde de Bethe (Masson, Paris, 

1983). 

M. Tomchenko, arXiv:cond-mat/1201.2185. 

C.E. Campbell, Phys. Lett. 44A, 471 (1973). 

E. Feenberg, Ann. Phys. 84, 128 (1974). 

LA. Vakarchuk and I.R. Yukhnovskii, Theor. Math. 

Phys. 40, 626 (1979). 

P.A. Whitlock, D.M. Ceperley, G.V. Chester, and M.H. 

Kalos, Phys. Rev. B 19, 5598 (1979). 

M.H. Kalos, M.A. Lee, P.A. Whitlock, and G.V. Chester, 



[20 
[21 

[22 

[23 

[24 

[25 
[26 
[27 

[28 
[29 
[30 
[31 



[32 
[33 

[34 



[35 



Phys. Rev. B 24, 115 (1981). 
E. Krotscheck, Phys. Rev. B 33, 
J.L. Epstein and E. Krotscheck, 
(1988). 



3158 (1986). 

Phys. Rev. B 37, 1666 



T. MacFarland, 
and M.H. Kalos 
L. Reatto, G.L. 
197, 189 (1994) 
I.O. Vakarchuk, 



S.A. Vitiello, L. Reatto, G.V. Chester, 
Phys. Rev. B 50, 13577 (1994). 
Masserini, and S.A. Vitiello, Physica B 



Rovenchak, J. 



V.V. Babin, and A. A. 
Phys. Stud. (Lviv) 4, 16 (2000). 
M. Tomchenko, Low Temp. Phys. 32, 38 (2006). 
R. Rota, J. Boronat, J. Low Temp. Phys. 166, 21 (2012). 
LA. Vakarchuk and I.R. Yukhnovskii, Theor. Math. 
Phys. 42, 73 (1980). 

MD. Tomchenko, Ukr. J. Phys. 50, 720 (2005). 

E. Zaremba and W. Kohn, Phys. Rev. B 15, 1769 (1977). 

K.V. Grigorishin, B.I. Lev, Ukr. J. Phys. 56, 1182 (2011). 

L.D. Landau and E.M. Lifshitz, Quantum Mechanics. 

Non-Relativistic Theory (Pergamon Press, New York, 

1980). 

S.A. Yushchenko, private communication. 

A.R. Jansen and R.A. Aziz, J. Chem. Phys. 107, 914 

(1997). 

M. Przybytek, W. Cencek, J. Komasa, G. Lach, B. 
Jeziorski, and K. Szalewicz, Phys. Rev. Lett. 104, 183003 
(2010); 108, 129902 (E) (2012). 

M.R. Andrews, D.M. Kurn, H.-J. Miesner, D.S. Durfee, 



15 



C.G. Townsend, S. Inouye, and W. Ketterle, Phys. Rev. 
Lett. 79, 553 (1997). 

[36] L.P. Pitaevskii, Sov. Phys. JETP 9, 830 (1959). 

[37] J.M. Ziman, Principles of the Theory of Solids (Cam- 
bridge Univ. Press, Cambridge, 1964), Chapter 1. 

[38] E.H. Lieb, W. Liniger, Phys. Rev. 130, 1605 (1963); E.H. 
Lieb, Phys. Rev. 130, 1616 (1963). 



[39] D. Ruelle, Statistical Mechanics. Rigorous Results (Ben- 
jamin Inc., New York, 1969). 

[40] I.M. Khalatnikov, An Introduction to the Theory of Su- 
perfluidity (Perseus, Cambridge, 2000), Chapter 1. 

[41] M. Tomchenko, arXiv:cond-mat/1201.2623. 

[42] M. Tomchenko, arXiv:cond-mat/1204.2149. 



